Phase estimation without a priori knowledge in the presence of loss 
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We find the optimal sclieme for quantum pliase estimation in tlie presence of loss when no a priori 
knowledge on the estimated phase is available. We prove analytically an explicit lower bound on 
estimation uncertainty, which shows that, as a function of number of probes, quantum precision 
enhancement amounts at most to a constant factor improvement over classical strategies. 
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I. INTRODUCTION 

Owing to highly promising predictions of the the- 
ory of precise quantum measurements and parameter 
estimation, as well as significant progress in quantum 
state engineering, the task of phase shift determin- 
ation has recently been readdressed both theoretic- 
ally and experimentally [1-11]. In classical systems 
the precision of the estimated phase scales with the 
amount of available resources as 1 / \/N, the so called 
Standard Quantum Limit (SQL) or more commonly 
the "shot noise" . Traditionally, iV denotes the num- 
ber of independent measuring probes, repetitions or 
copies of a system. The potential precision boost 
offered by quantum mechanics stems from the pos- 
sibility of preparing TV copies of a system in a highly 
entangled state, particularly sensitive to the vari- 
ations of the estimated parameter [1-3]. In ideal 
scenarios, these states yield phase estimation preci- 
sion which scales as 1/N and is referred to as the 
Heisenberg Limit (HL). 

Environmentally induced decoherence, however, 
significantly affects the performance of entanglement 
based quantum strategies [12-21] with photon loss 
being its most relevant source in optical implementa- 
tions. The need to balance the phase sensitivity and 
robustness against losses results in states perform- 
ing better than SQL yet falling short of HL [18, 19]. 
Other approaches, trying to mimic the quantum en- 
hanced strategies using multiple-pass technique [8] 
are even more susceptible to losses and cannot com- 
pete with the optimally designed entangled states 
[22]. Despite the quantitative improvement of pre- 
cision offered by quantum states in the presence of 
loss, it has remained an unsolved problem whether in 
the asymptotic regime N oo quantum states offer 
better than SQL scaling, i.e., c/N" with a > 1/2. 

In this paper we solve the problem of optimal 
phase estimation in the presence of loss with no a 
priori knowledge, and prove analytically that even 
for arbitrarily small loss, quantum enhancement does 
not offer better than c/\/N scaling for N — > oo, and 
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Figure 1. Phase estimation setup. Channel a acquires a 
phase (fi relative to channel b. Losses are modeled by two 
beam splitters with power transmissions rja, t?6. 



the only gain over classical strategies is a smaller 
multiplicative constant c. It should be emphasized 
that the proof contains the most general description 
of a quantum measurement, hence its conclusions 
are valid also for adaptive schemes (see Appendix 
C), which are especially interesting from a practical 
point of view [23, 24]. 



II. MODEL 

Two approaches to phase estimation are typically 
pursued. In the first, local approach, a measurement 
scheme is devised, which offers the highest sensitiv- 
ity to phase deviations from an a priori known value, 
If = (fQ. This is achieved by finding a strategy 
that maximizes the quantum Fisher information, 
Fq, which defines the lower bound on the precision 
of the estimated phase through Sip > 1/ y/Fq [25- 
28]. The optimal states have been found both for 
lossless [1, 2] (the so called NOON states) and more 
realistic lossy scenarios [18, 19]. 

The second approach, which we will pursue in this 
paper and refer to as the global approach, assumes 
no a priori knowledge about the phase, so that if is 
equiprobably distributed over the [0, 27r) region. We 
consider a general pure N photon two-mode state 
[29] 
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which is fed into an interferometer with a relative 
phase delay (see Fig. 1). Apart from acquiring the 
phase via the unitary U^^ = ^^^'^'^^'^^ the state ex- 
periences losses modeled by two beam splitters with 
power transmissions ryo and rib [30] . The output 
state then takes the form Pout(<p) = U^poutU^, where 

pout = E E \^'''''){^'-'% (2) 

with subnormalized conditional states corresponding 
to la and lb photons lost in arms a and b respectively 

N-h 

1^'-'")= J2anPl:'''\n-la,N-n-lb) (3) 

n=la 

where 

= ^jBiMB^^-^inb), Bnv) = 0(1-^)''?" 

(4) 

Keeping the reasoning most general, the information 
about ip is extracted via a measurement on pouti'p) 
described by a Positive Operator Valued Measure 
(POVM), {Mr}, Y^r^r = 1- The outcome r is 
observed with probability p {r\ip) = Tr {pout{v)Mr}, 
and the estimated phase inferred from it is defined 
by an estimator (p (r) . Optimization procedure with 
respect to a given cost function C{ip, (p) amounts to 
finding the state the measurement {My.}, and 
the estimator (p (r) that minimize the cost function 
averaged over a flat a priori phase distribution 

Let ^) = = Er=-oo c„e'"('^-'^), be an 

arbitrary real symmetric cost function {cn = C—n < 
for n ^ 0) respecting the cyclic nature of ip [26, 31]. 

III. OPTIMIZATION 

Thanks to the flat a priori phase distribution, the 
problem enjoys a symmetry with respect to an ar- 
bitrary phase shift U^. The search for the optimal 
measurement strategy may be restricted to the class 
of covariant POVM {M^} [26, 31, 32] parameterized 
by a continuous parameter cp: M(p — UipSU^, where 
S is a positive semi-definite operator satisfying the 
POVM completeness constraint J ^U^pSU^ = 1. 
With the above substitution, the average cost func- 
tion simplifies to 

(C)= I gTr{pout(v')S}C(^) (6) 



and (C) has to be minimized only over the choice of 
the input state \4>)in and the seed operator S. 

In order to find the optimal S, one can re- 
write Eq. (2) in the form Pout = ®N'=oPout^ with 

Po^u't = E^Zo' \^'-''-'''-''H^'-''-'''-'^l which 

reveals the block structure with respect to the total 
number of surviving photons A^'. Therefore, without 
loss of generality, we may impose an analogous block 
structure on the seed operator S = ® j^/_o ■ 
Physically, such a block structure implies that a non- 
demolition photon number measurement had been 
performed at the output, before any further phase 
measurements have taken place. Following the reas- 
oning presented in [26, 31] it can be shown that 
without loosing optimality, the input state para- 
meters an can be chosen real, in which case the 
optimal seed operator S^^. = |ejv') (cat'I, where 

\eN') = J2n=o l'^' ^' ~ Appendix A). 

In what follows we choose the cost function 
C{ip-^) = 4sin2^ (co = 2,ci = c_i = -1) 

and denote its average by 5^<p, as it is the simplest 
cost function approximating the variance for narrow 
distributions [3]. 

Performing the integration in Eq. (6) the average 
cost function reads: 

5^ = 2- a+Aa, (7) 

where non-zero elements of the matrix A read: 

n,N—n 

Hence, the minimal cost equals 5'^p = 2 — Amax) 
where Amax is the maximal eigenvalue of the mat- 
rix A, and the corresponding eigenvector provides 
the optimal input state parameters a. 



A. Numerical solution 

Numerical results of the above eigenvalue prob- 
lem are presented in Fig. 2. Black lines depict phase 
estimation uncertainty 6^p of the optimal quantum 
strategy plotted as a function of N for r]a = Vb G 
{0.6,0.8,1}. In the absence of loss the optimal 
quantum curve tends to the Heisenberg scaling, 
whereas, when losses are present, it flattens signi- 
flcantly with increasing N. The inset depicts the 
form of the optimal state. With increasing degree 
of loss the distribution of a„ for the optimal state 
becomes more peaked as compared with the lossless 
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Figure 2. Log-log plot of optimal phase estimation uncer- 
tainty as a function of number of photons used for three 
different levels of loss (equal in both arms): rj = 1 (solid), 
rj — 0.8 (dashed), rj — 0.6 (dotted). White area in the 
middle of the picture corresponds to 1/A'' < 5^3 < 1/y/N. 
Gray lines represent asymptotic bounds given by Eq. (12) 
for r] = 0.8, rj = 0.6. The inset depicts the structure of 
the optimal states for the three levels of loss for A'^ — 100. 
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{n+l)T 
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[3] . This behavior can 

be intuitively understood in a similar fashion as in 
the local approach [18, 19], where the NOON states 
with only two non-zero coefficients ao^a^, are the 
most sensitive to the phase shift but extremely vul- 
nerable to loss. In the presence of loss, larger weights 
need to be ascribed to intermediate coefficients, in 
order to preserve quantum superposition even after 
some photons are lost. The same effect of increasing 
weights of intermediate coefficients at the expense of 
marginal ones is also present in the global approach. 



v^Av < Amax- Let a be the eigenvector correspond- 
ing to Ainax^ cc^Aa = Amax- The fact that all matrix 
elements of A are non-negative, implies V„ck„ > 0. 

Let us now define a matrix A', such that all 
nonzero entries of A are replaced by the maximum 
matrix element = max„ {A„.„_i} = Ajv,Ar-i. 
Since a„ > and A' „, > An m > we can write: 



An 



(9) 



where X'^i^^ is the maximal eigenvalue of A'. X'^^^^^ 
can be found analytically by noting the following re- 
currence relation for the characteristic polynomial of 
A': detA„+i — — AdetA„ — A''^^ det A„^i, where 
A = A' — Al, while A„ are (n -I- 1) x (n -I- 1) 
submatrices of A. The solution of the recurrence 
relation reads det (A) = Z?Ar_|_i(— A, yl''^^), where 



Dn{x,a'^) = a' 



, sin[(n+i)arccos(^)] 



is the Dickson 



sin[arccos(^)] 
polynomial [33] of the nth order. The largest ei- 
genvalue corresponds to the largest root of det (A), 

A' „^ — 2A'^cos 



{N+2) 

We can finally write explicitly the lower bound on 
the variance: 



S^ip > 2 



^-cos(j^)f:^B^i,)Br\v) 



1=0 



(10) 

Expanding the above formula in the limit N ^ oo 
we get: 



5V > 



1-^ 
4t]N 



O 



1 

iV2 



(11) 



B. Asymptotic bounds 

We now move on to present the main result of the 
paper. Numerical results presented above and the 
ones obtained within the local approach [18, 19] in- 
dicate that in the presence of loss, phase estimation 
uncertainty Sep departs from the HL and asymptotic- 
ally approaches c/^/N. Until now, however, an ana- 
lytical proof of the above conjecture was missing. 

Let us first derive an upper bound on the maximal 
eigenvalue Amax of matrix A in Eq. (7). Without 
loss of generality, we assume that r]a < i]h. Clearly, 
setting r]h — I can only improve our estimation — 
hence Amax increases. For rja = rj < 1, rji, = 
1 the nonzero matrix elements read: A„.„_i = 

Recall that for an arbitrary normalized vector v, 



which proves that for rj < 1, S(p scales as c/^/N. 

A tighter bound can be analogously derived for the 
case ria = rib = rj, by noting that max„ {A„^„_i} = 
A|-jv-| |"iv-|_.^. In the limit iV — ^ oo we get: 



C. Optimal classical strategy 

For the sake of comparison, we also derive the op- 
timal classical phase estimation strategy, in which a 
coherent state with mean photon number N is sent to 
an initial beam splitter of transmissivity Tin, whose 
output feeds paths a and b of the interferometer. 
We assume no additional external phase reference, 
hence the state is effectively a mixture of terms with 
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a different total photon number. The optimal seed 
POVMis e?^,^oS^p; yielding: 

2B[Ni^an^]B[N7^,(l-n^)\ _ ^^^^ 

Ny^riaTinVb (1 - Ti„) 

where B{x) = e~'^ fjv^ is the Bell polyno- 

mial of order 1/2. For strong beams {N — >• oo) up to 

the first order in 1/N, 5^ « + ^i^^) 

and is minimized for the choice Tin = 1/(1 + -s/Va/Vb) 



1 1 



^ ' (14) 



which is exactly the same formula as for the optimal 
classical strategy in the local approach [19] . 



amounts just to a smaller multiplicative constant c 
in the scaling law c/^/N. Comparing Eq. (14) (with 
Va = Vb = 1) with the bound given in Eq. (11) 
we may conclude that asymptotically quantum en- 
hanced protocols provide at most a factor of 




decrease in the uncertainty of estimation. In the 
case ria = rjb = r], using a tighter bound (12) the 
above factor reads 1/^1 — rj. We conjecture that the 
fact that losses necessarily turn HL into c/ a/JV is a 
general feature of all quantum estimation problems, 
such as estimation of direction, Cartesian frames etc. 
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IV. CONCLUSIONS 

Results presented in the paper indicate that, while 
quantum enhanced protocols provide quantitative 
boost in the estimation precision, the presence of loss 
unavoidably causes the precision scaling to become 
classical in the limit of large number of resources N. 
The asymptotic gain of quantum enhanced protocols 
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Appendix A: Optimal measurement 



Substituting the output state p^ut — ®jv'=o Pout ^^'^ ^^^'^ operator S = 0]^/^o ^° (^)' ^® S^t 
an explicit formula for the average cost function: 



N N-N' N'+la. 

= E E E 

JV'=0 la=0 n.m- 



(Al) 



where k = N - N' - la, C„m = / ^C(<^)e*("-™)^ and S^',„, = (n',iV' - n'jS^'jm',^' - m'). The 

completeness constraint / ^U^EU^ = t implies that = 1. Therefore, if restricted to m = n terms, the 

sum (Al) reduces to a constant cq = Coo- Changing the summation order we can rewrite Eq. (Al) as 



N min{n,7n) A^ — max(n,7n) 



n,m— la—0 



h=0 



'n — la,Tn — la' 



(A2) 
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Now, as for all n m cost coefficients Cnm < 0, we get the following lower bound on the average cost 

TV inin(7i,m) N — ma^x(n.m) 

{C)-Co> E E Cn^/S^T''^/?^''" Kl l^ml 

n,m=0 ;„=0 lh=0 

N min(n,m) AT— max(n,m) 

^ E E E Cr^mPy'Pl^'''\a:\\am\. (A4) 

ra,m=0 la=0 lb=0 

The first inequality is saturated by choosing input state's and seed operator's coefficients to be real. The 
second inequality follows from '^^',m' — \J'^m,',m''^n'',n' ~ which is a consequence of positive semi- 
defniteness of and the completeness constraint. Both inequalities are saturated for S^,. = |ejv') (ejv'l) 
where |ejv') = X^^=o l^^i ^' ~ This proves the optimality of the measurement considered in the paper. 



(A3) 



Appendix B: Distinguishability of photons 

If photons traveling through the interferometer are distinguishable, e.g. they are prepared in different time 
bins, the dimension of the Hilbert space needed to describe the state of N photons is 2^, as opposed to N+1 
for the indistinguishable case. In fact, the indistinguishable case may be considered as a restriction of the 
former space to its fully symmetric subspace. We prove below that considering distinguishable photons is of 
no use, since the optimality can always be attained within the class of states belonging to the fully symmetric 
(bosonic) subspace. Let 

iV'iv) = E I") ' (^^) 

n=0'' 

be a general state of N distinguishable photons traveling through the interferometer, where the sum runs 
over all A''-bit sequences n, with |n) = jm) (g) . . . In^v), where \ni) = |1) (|0)) denotes a photon in the ith time 
bin, propagating in the a{b) arm of the interferometer respectively. 

Taking loss into account, we additionally need to track the time slots in which photons were lost. We 
define a binary string la = la,il'a,2 ■ ■ ■ la,N with Is representing the time bins in which photon was lost in 
arm a and similarly h for the arm b. The general seed operator has a block diagonal structure with respect 

to different patterns of surviving photons: = ®jv'=o" " ' where Is in the binary string TV denote the 
time bins in which photons were successfully transmitted. Formally, using bitwise subtraction, we can write 

N' = 1 — la — lb- Written in a basis = X]7i',m'=o»' which n' stands for a string with 

A^' bits placed at positions corresponding to Is in AT' with complementary positions left empty (neither 
nor 1). In order to simplify the notation, for any binary sequence x, we denote by a; = \x\ the number of 
Is in the sequence. Moreover, we use a notation a; \ y for a binary string x with empty entries at positions 
corresponding to Is in y. 

Adapting Eq. (A2) to the distinguishable photon case, we get: 

1 min{Ti.TTi-) 1 — max(Ti,m) 

{C)-co= Yl E E ^-7i-S^^'^<ar„S^\£:;:)U\(..+..) (B2) 

n,m=0 la=0 lb=0 



where min, max should be understood as bitwise operations, 7^°''*' = y (1 — rjaY^ria "(1 — ribY'"!]^ " and 
for simplicity we have put = 0^, 1 = 1^. 
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We now split the sums over li into sum over (number of Is in Zj) and the sum over permutation of Is 
within li. We proceed analogously for summations over n (m) obtaining 

AT min(ri,m) A^ — max(n,m) 

{C)-co= E E C'„™7i-S^^''' (B3) 



1 min(Ti,m) 1 — max(Ti,m.) 

"Ti\(J<,+Z6),r7i\(J„+ii, 



E E E E 



n=0 Tn.=0 ia.=0 i(,=0 

\nl=nlml = m llal=la |ibl=ib 



In order to proceed further let us for the moment specialize to lossless case r/o = % = 1, where the above 
formula simplifies to: 



Af 11 

(C) - CO = 53 C„„ ^ ^ a^ar^S^ (B4) 

n,m=0 n=0 m=0 

n^im |»i|=n |m.| = m 

S needs to be a positive semi-definite operator, and by completeness constraint Em^n = <^m,,n7 whenever 
n = m. Since diagonal blocks of S (corresponding to n = to) arc proportional to identity, it implies that 
none of the off-diagonal blocks of S (corresponding to n ^ m) can have a singular value larger than 1. This 
can be proven as follows. Let us assume that for certain block (m, n) (n ^ m), the largest singular value 
A > 1, and let \vm), \wn} be the normalized left and right singular vectors corresponding to singular value 
I'Wml = m, \wn\ = n. Defining \z) = \vm) — \wn), we calculate 

{z\E\z) = (-u„|S|t;„) + {w^\E\w^} ~ 23?(t;„|S|u;„) = 2(1 - A) < 0, (B5) 

which contradicts the positivity semi-definiteness of S. Because all singular values of any (n,m) block of S are 

smaller than one, the following inequality holds: X]n=o S to=o o:^o:mE.n,m < a*^oim, = /X^ti^o 

|n|=,i |m|=m Y \-n.\=n 

This leads to a bound on the cost function in the lossless case 

N 

(C) - Co > ^ CnmOLnO*^, (B6) 
n,m— 

proving that one can achieve optimality restricting oneself to indistinguishable photons. 

Returning to Eq. (B4), we see that we can apply a similar argumentation making use of positive semi- 
definiteness of S^,'^''^ = X^™^^o ""' E^^^^^i^i^^ n\(j„+ii,) operator. We notice that the complete- 

ness constraint again implies a block structure of S^'"''''^ with respect to to = |m|,n = |n|, with diagonal 
elements of diagonal blocks (n, n) being now Y^^i^=o J2\^=o ^ ~ {?) i^f")- "^^^^ implies that the maximum 

singular value of any (m,n) block of is constrained by (™"(^''")) -ma.K{n,m)y ^ ^.^^^j^^ ^^^^^jj^ 

the following bound: 

min(n,m) A^ — max(n.rn) / . / \\ /at / \\ 



ra,m=0 la=0 lb=0 

n^ra 



N min(rt,m) AT— max(n,m) 

E E E 

n,m=0 ;<,=0 ^6=0 



(B7) 



E E E '^.™*'"^^N/(;:)(i')r,:")f"r')KiK 



.RecalUng that = (;") (^i^")')'" ''''' evident that the above equation is identical to Eq. (A4) 

obtained for the indistinguishable case, this completes the proof that the optimal estimation is indeed 
achievable using indistinguishable photons. 
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Appendix C: Adaptive measurement schemes 

Let us describe a general structure of adaptive measurement schemes performed on N subsystems. Let 
{n^^^} be a POVM performed on the first copy. Depending on the measurement result ii a POVM {11,-^'' (*i)} 

is performed on the second copy. In general, a POVM performed on the fc-th copy {n-^^(ii, . . . ,ik-i)} depends 
on all previous measurement results. The adaptive measurement mathematically corresponds to POVM: 

ni = n,„...,,„ =n«®---®n(;;)(ii,...,i;v-i), (ci) 

where Xlj can be treated as a single global POVM with measurement results indexed by i. This shows 
that, for distinguishable subsystems, optimization of estimation strategy over global POVMs covers also the 
case of adaptive measurements. Moreover, we have proved earlier in Appendix B, that the optimal phase 
estimation can be realized using indistinguishable subsystems. Therefore, the bounds derived in the paper, 
which assume a global POVM on indistinguishable photons, indeed hold also for all adaptive measurement 
strategies. 
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